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Localization of electromagnetic waves in a two dimensional random medium 
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Motivated by previous investigations on the radiative effects of the electric dipoles embedded in 
structured cavities, localization of electromagnetic waves in two dimensions is studied ab initio for 
a system consisting of many randomly distributed two dimensional dipoles. A set of self-consistent 
equations, incorporating all orders of multiple scattering of the electromagnetic waves, is derived 
from first principles and then solved numerically for the total electromagnetic field. The results show 
that spatially localized electromagnetic waves are possible in such a simple but realistic disordered 
system. When localization occurs, a coherent behavior appears and is revealed as a unique property 
differentiating localization from either the residual absorption or the attenuation effects. 
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When propagating through a medium consisting of 
many scatterers, waves will be scattered by each scat- 
terer. The scattered waves will be again scattered by 
other scatterers. Such a process will be repeated to es- 
tablish an infinite recursive pattern of multiple scatter- 
ing. As a result, the wave propagation may be signif- 
icantly altered!! It is now well-known that the multi- 
ple scattering of waves is responsible for many fascinat- 
ing phenomena, ranging from phenomena of macroscopic 
scales such as twinkling lights in the evening sky, modula- 
tion of ambient noise in the oceans, and electromagnetic 
scintillation of turbulence in the atmosphere, to phenom- 
ena oLmicroscopic or mesoscopic scales including random 
lasersB and electronic resistivity in disordered solidsa. It 
has also been proposed that under certain conditions, the 
multiple scattering can lead to the unusual phenomenon 
of wave localization, a concept originally introduced to 
describe disorder induced metal-insulator transitions in 
electronic systems!! 
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calization phenomenon has been characterized by two 
levels. One is the weak localization associated with the 
enhanced backscattering. That is, waves which propa- 
gate in the two opposite directions along a loop will ob- 
tain the same phase and interfere constructively at the 
emission site, thus enhancing the backscattering. The 
second is the strong localization, without confusion often 
just termed as localization, in which a significant inhi- 
bition of transmission appears and the energy is mostly 
confined spatially in the vicinity of the emission site. 

While the weak localization, regarded as a precursor to 
the strong localization, has beenr-well studied both the- 
oretically (e. g. the monographBEJ) and experimentally 
(e. g.Ej), observation of strong localization of classical 
waves for, higher than one dimension remains a subject 
of debatcotiJEj, primarily because a suitable system is 
hard to find and the observation is often obscured by 
such effect &-as the residual absorption!^ and scattering 
attenuatiorcj. 



In this paper, we wish to present a simple but realistic 
system to study the phenomenon of strong localization 
in two dimensions. The system we consider stems from 
the previous research on enhancement and inhibition of 
electromagnetic radiation in structured media. For ex- 
ample, KuhrcJ used the model to describe the effect of 
a metallic mirror on the radiatipa from a nearby excited 
molecule, whereas Chance et al.Ej used to explain the ex- 
perimental d§ta of Drexhage. Later, Erdogan et alBa and 
Wang et alE3 employed the model to analyze the effects 
of the cylindrically periodic structure on the radiation 
of an enclosed two-dimensional (2D) dipole. Inspired by 
the work of Erdogan et al.ES, we attempt to construct 
a system consisting of many 2D dipoles. The propaga- 
tion of electromagnetic (EM) waves in such a system is 
formulated rigorously in terms of a set of coupled equa- 
tions, and then is solved numerically. We show that a 
strong localization of EM waves is possible in this.-s.ys- 
tem. In line with the work on acoustic localizations, a 
phase diagram method is used to describe localization of 
EM waves. Since the system considered here results from 
the practical application of light emission in structures, 
an experimental verification-may be readily realized. 

Following Erdogan et al.c3, we consider 2D dipoles as 
an ensemble of harmonically bound charge elements. In 
this way, each 2D dipole is regarded as a single dipole 
line, characterized by the charge and dipole moment per 
unit length. Assume that N parallel dipole lines, aligned 
along the z-axis, are embedded in a uniform dielectric 
medium and randomly located at fi (i = 1, 2, . . . , N). The 
averaged distance between dipoles is d. A stimulating 
dipole line source is located at r s , transmitting a contin- 
uous wave of angular frequency u). By the geometrical 
symmetry of the system, we only need to consider the z 
component of the electrical waves. 

Upon stimulation, each dipole will radiate EM waves. 
The radiated waves will then repeatedly interact with the 
dipoles, forming a process of multiple scattering. The 
equation of motion for the z-th dipole is 
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where cjqj is the resonance (natural) frequency, pi, qt and 
rrii the dipole moment, charge and effective mass per unit 
length of the i-th dipole respectively. E z (fi) is the total 
electrical field acting on dipole pi, which includes the ra- 
diated field from other dipoles and also the directly field 
from the source. The factor bo,i denotes the damping and 
is determined by energy conservation. Without dissipa- 
tion, &o,i can be determined from the balance between 
the radiative and vibrational energies, and is given asc3 
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with e being the constant permittivity and c the speed of 
light in the medium separately. .— . 

Equation (|l|) is similar to Eq. (1) ir£3. The difference 
lies in the driving field on the right hand side of the 
equation. Inc3, E z is the reflected field at the dipole due 
to the presence of reflecting surrounding structures, while 
in the present case the field is from the stimulating source 
and the radiation from all other dipoles. 

The transmitted electrical field from the continuous 
line source is determined by the Maxwell equationsBj 
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where u> is the radiation frequency, and po is the source 
strength and is set to be unit. The solution of Eq. (g) is 
clearly 



G (r~ f* s ) = {p, bJ 2 po)iitH^\k\r - r a \)e~ 



(4) 



with k — lu/c, and Hq 1 ^ being the zero-th order Hankcl 
function of the first kind. 

Similarly, the radiated field from the i-th dipole is 
given by 



V 2 - 



d 2 



2 dt 2 



Gi(f- n) = » SL Pi 5( 2 \r- n). (5) 



The field arriving at the i-th dipole is composed of the 
direct field from the source and the radiation from all 
other dipoles, and thus is given as 



E z (fi) = G (n - r s ) + 



N 

E 



Substituting Eqs. (|), (|), and 
writing p, = p t e~ 



Gift-?,). (6) 
* into Eq. (pi), and 
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This set of linear equations can be solved numerically for 
Pi. After pi are obtained, the total field at any space 
point can be readily calculated from 



N 
3 = 1 



Pi)- 
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In the standard approach to wave localization, waves 
are said to be localized when the square modulus of the 
field |i?(r)| 2 , representing the wave energy, is spatially 
localized after the trivial cylindrically spreading effect is 
eliminated. Obviously, this is equivalent to say that the 
further away is the dipole from the source, the smaller its 
oscillation amplitude, expected to follow an exponentially 
decreasing pattern. 

To this end, it is instructive to point out that an al- 
ternative two dimensional dipolejnodel was deviced pre- 
viously by Rusek and Orlowskilla. The authors derived 
a set of linear algebraic equations, which is similar in 
form to the above Eq. (]?]). However, there are some fun- 
damental differences between the two models. InE_3, the 
interaction between dipoles and the external field is de- 
rived by the energy conservation, while in the present 
case the coupling is determined without ambiguity by 
the Newton's second law. The former leads to ar. 
determined phase factor. According to, e. g. Refs.t^l 
the energy conservation can only give the radiation fac- 
tor in Eq. (|^). We would also like to point out that 
the set of couple equations in Eq. (Q) is similar in spirit 
to the tighUbkiding model used to study the electronic 
localizationclE3. 

There are several ways to introduce randomness to 
Eq. (0). For example, the disorder may be introduced by 
randomly varying such properties of individual dipoles as 
the charge, the mass or the two combined. This is the 
most common way that the disorder is introduced into 
the tight-binding model for electronic wavesEa. In the 
present study, the disorder is brought in by the random 
distribution of the dipoles. 

Before moving to solve Eq. (|t]) for the phenomenon of 
localization of EM waves, we discuss a general property 
of wave localization. The energy flow of EM waves is 
J ~ E x H . By invoking the Maxwell equations to relate 
the electrical and magnetic fields, we can derive that the 
time averaged energy flow is 



< J >t= 



dtJ - |£rV6> 



(9) 



where the electrical field is written as E = eE\E\e ie , with 
ce denoting the direction, |£?| and 9 being the amplitude 
and the phase respectively. It is clear from Eq. (||) that 
when 9 is constant, at least by spatial domains, while 
\E\ ^ 0, the flow would come to a stop and the energy 
will be localized or stored in the space. In the localized 
state, a source can no longer radiate energies. Alter- 
natively, we can write the oscillation of the dipoles as 
Pi = \Pi\e l8i ■ By studying the square modulus of pi in 
the form of | f* j — r s \ \pi \ 2 , and its phase 9i , we can also 
investigate the localization of EM waves. Note here that 
the factor \ fi — f s \ is to eliminate the cylindrical spread- 
ing effect in 2D as can be seen from the expansion of 

the Hankel function \Hq 1 \x)\ 2 ~ -. That the phase 9 
is constant implies that a coherence behavior appears in 
the system. In other words, the localized state is a phase- 
coherent state. This picture of localization also holds for 
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acoustic and electronic waves. It is a unique feature of 
wave localization. 

For simplicity yet without losing generality, assume 
that all the dipoles are identical and they are randomly 
distributed within a square area. The source is located at 
the center (set to be the origin) of this area. For conve- 
nience, we make Eq. (^) non-dimensional by scaling the 
frequency by the resonance frequency of a single dipole 
luq. This will lead to two independent non-dimensional 

parameters b — and b' = ^ (f^)' Both parame- 
ters may be adjusted in the experiment. For example, the 
factor bo can be modified by coating layered structures 
around the dipolescj. Then Eq. (uu becomes simply 
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with / = Eq. (|T(]) can be solved numerically for 
and then the total field is obtained through Eqs. (|), (P) 
and (ph. In the calculation, we scale all lengths by the 
averaged distance between dipoles d. In this way, the fre- 
quency to always enters as kd and the natural frequency 
ujo as kod. We find that all the results shown below are 
only dependent on parameters b, bo/tuo, and the ratio 
lo/ujo or equivalently fc/fco = (kd) / (kod) . Such a simple 
scaling property may facilitate designing experiments. 

We have first computed the transmitted intensity av- 
eraged over the random distribution of the dipoles as a 
function of non-dimensional frequency kd. The results in- 
dicate that the natural frequency of the dipoles is altered 
by the multiple scattering of EM waves, and is shifted 
towards the lower frequency end. We also find that the 
transmission is significantly suppressed in a range of fre- 
quencies slightly above the natural frequency of a single 
dipole, suggesting the strong localization in this range 
of frequencies. This helps us to search for the frequen- 
cies where the strong localization is most prominent. We 
show an example below. 

By setting b = bo/^o — 10~ 3 and kod — 1, we find a 
strong localization of EM waves at, for instance, lo/loo = 
1.008. In this case, the natural frequency of the dipoles 
is shifted to around 0.9<x>o- A typical picture of wave 
localization is shown in the bottom part of Fig. [l] for an 
arbitrary random distribution of 400 dipoles. To describe 
the phase behavior of the system, we assign a unit phase 
vector, Ui — cos^e^ + sinf^ej, to the oscillation phase 
6i of the dipoles. Here e x and e y are unit vectors in the 
x and y directions respectively. These phase vectors are 
represented by a phase diagram in the x — y plane with 
the phase vector Hi being located at the dipole to which 
the phase 0i is associated. The phase behavior of the 
system is depicted by the top portion of Fig. [l| 

Here it is clearly shown that the field energy is strongly 
localized near the transmitting source, and, as expected, 
the energy decreases almost exponentially along the ra- 
dial direction. Meanwhile, the system reveals an in-phase 
phenomenon: nearly all the phase vectors of the dipoles 
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FIG. 1: Top: The phase diagram for the two dimensional 
phase vectors denned in the text; the phase of the source is 
assumed to be zero. Bottom: The spatial distribution of en- 
ergy (~ l^| 2 )- Here the geometrical factor has been removed. 



point to the same direction, exactly opposite to the phase 
vector of the source, i. e. the dipoles tend to oscillate out 
of phase with the source. It is easy to see that such a 
coherent behavior effectively prevents waves from prop- 
agation. The picture represented by Fig. [I] fully com- 
plies with the general description of wave localization 
stated above. The energy localization in Fig. |l| is also 
in qualitative agreement with that observed in the mi- 
crowave localizatonE3E£ -and the acoustic localization in 
water with air-cylindersca. 

We note from Fig. [l] that near the sample boundary, 
the phase vectors start to point to different directions. 
This is because the numerical simulation is carried out 
for a finite sample size. For a finite system, the energy 
can leak out at the boundary, resulting in disappearance 
of the phase coherence^]. When enlarging the sample 
size by adding more dipoles while keeping the averaged 
distance between dipoles fixed, we observe that the area 
showing the perfect phase coherence will increase accord- 
ingly. Another factor affecting the phase coherence be- 
havior, and subsequently the wave localization, is the 
damping factor. When the absorption is added, the in- 
phase behavior will be degraded, and the waves will be- 
come delocalized gradually, in agreement with the previ- 
ous simulation on acoustic localization in bubbly watero. 
Though degrading, we find that the localization can sus- 



4 




r/d 

FIG. 2: Energy versus the distance away from the source. 
Again we removed the cylindrical spreading factor. 

tain a substantial variation in the damping factor, mak- 
ing the system a good candidate for observing the local- 
ization phenomenon. 

To find the localization length in Fig. [j], we plot the 



total energy in all directions as a function of the distance 
from the source. The results are presented by Fig. ||. 
Here, the numerical data are shown by the black squares, 
and the result fitted from the least squares method is 
shown by the solid line. It shows that after removing 
the spreading factor, the data can be fitted by e - '^. 
From the slop of the solid line, the localization length £ 
is estimated as around 2.02d, which is in the vicinity of 
the localization length 1.6d estimated experimentally for 
microwave localization in 2D dielectric latticesEl 

In summary, we have presented a model system to 
study the localization of EM waves in 2D random me- 
dia. The results show that spatially localized EM waves 
can be realized in this simple but realistic disordered sys- 
tem. When localization occurs, a coherent behavior ap- 
pears as a unique feature separating the localization from 
other effects such as the residual absorption or extinction 
effects. Such a coherence phenomenon for waves in ran- 
dom media could also be relevant to the random laser 
actionnl. 

This work is supported by the Bai Yu Lan Fund of 
Shanghai and the NSF of China (20074007, 90103034). 



* Corresponding author: Department of Physics, Na- 
tional Central Unive rsity, Taiwan; Electronic address: 
2hen@phy.ncu.edu. tw 

1 A. Ishimaru, Wave propagation and scattering in random 
media, (Academic Press, New York, 1978). 

2 N. M. Lawandy, R. M. Balachandran, A. S. L. Gomes, and 
E. Sauvain, Nature 368, 436 (1994). 

3 S. Datta, Electronic transport in mesoscopic systems, 
(Cambridge, New York, 1995). 

4 P. W. Anderson, Phys. Rev. 109, 1492 (1958). 

5 S. John, Phys. Rev. Lett. 53, 2169 (1984). 

6 C. H. Hodges and J. Woodhouse, J. Acoust. Soc. Am. 74, 
894 (1983). 

7 T. R. Kirkpatrick, Phys. Rev. B 31, 5746 (1985). 

8 S. He and J. D. Maynard, Phys. Rev. Lett. 57, 3171 (1986). 
" C. A. Condat, J. Acoust. Soc. Am. 83, 441 (1988). 

D. Sornette and B. Souillard, Europhys. Letts. 7, 269 



10 



11 A. Z. Genack and N. Garcia, Phys. Rev. Lett. 66, 2064 
(1991). 

12 R. Dalichaouch, J. P. Armstrong, S. Schultz, P. M. Platz- 
man and S. L. McCall, Nature 354, 53 (1991). 

13 S. L. McCall, P. M. Platzman, R. Dalichaouch, D. Smith, 
and S. Schultz, Phys. Rev. Lett. 67, 2017 (1991). 

14 A. Lagendijk and B. A. van Tiggelen, Phys. Rep. 270, 143 
(1996). 

15 M. Rusek and A. Orlowski, Phys. Rev. E 51, R2763 (1995). 

16 M. M. Sigalas, C. M. Soukoulis, C.-T. Chan, and D. 
Turner, Phys. Rev. B 53, 8340 (1996). 

17 A. A. Asatryan, et al. Phys. Rev. B 57, 13535 (1998). 



23 



P. Sheng, Introduction to wave scattering, localization 
and mesoscopic phenomena, (Academic Press, New York, 
1995). 

Y. Kuga, A. Ishimaru, J. Opt. Soc. Am. Al, 831 (1984); 
M. van Albada, A. Lagendijk, Phys. Rev. Lett. 55, 2692 
(1985); P. E. Wolf, G. Maret, Phys. Rev. Lett. 55, 2696 
(1985); A. Tourin, et al., Phys. Rev. Lett. 79, 3637 (1997); 
M. Torres, J. P. Adrados, F. R. Montero de Espinosa, Na- 
ture 398, 114 (1999). 

D. S. Wiersma, P. Bartolini, A. Lagendijk, and R. Roghini, 
Nature 390, 671 (1997). 

A. A. Chabanov, M. Stoytchev, and A. Z. Genack, Nature 
404, 850 (2000). 

F. Scheffold, R. Lenke, R. Tweer, and G. Maret, Nature 
398, 206 (1999); D. Wiersma, et al. Nature, 398, 207 
(1999). 



Y.-Y. Chen a nd Z. Ye, to appear in Phys. Rev. E; |cond- 
Lt/0109494 



z4 H. Kuhn, J. Chem. Phys. 53, 101 (1970). 

25 R. R. Chance, A. Prock, and R. Silbey, in Advances in 
Chemical Physics, I. Prigogine and S. A. Rice, eds. (Wiley, 
New York, 1978), Vol. 37, pp. 1-65. 

26 T. Erdogan, K. G. Sullivan, and D. G. Hall, J. Opt. Soc. 
Am. B 10, 391 (1993). 

27 C.-C. Wang and Z. Ye, Phys. Stat. Sol. (a) 174, 527 (1999). 

28 E. Hoskinson, and Z. Ye, Phys. Rev. Lett. 83, 2734 (1999). 

29 R. Zallen, The physics of amorphous solids, (John Wiley 
and Son, New York, 1983). 

30 K.-X. Wang and Z. Ye, Phys. Rev. E 64, 056607 (2001). 



